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I . Introduction 

Airplanes flying at supersonic speed produce "sonic booms" on the 
ground. Sonic boom reflected on the ground is felt by people as a "start- 
ling noise," and the present levels and shapes of the overpressures are 
highly objectionable. The sonic boom is unavoidable for vehicles having 
lift and cruising in horizontal flight. Therefore, the direction to follow 
in order to attempt to minimize the unacceptability on the part of the 
people is to decrease the main source of annoyance which is largely con- 
nected with the presence of shock waves. 

The results of an analytical and experimental investigation are pre- 
sented and discussed herein, and have been performed in order to determine 
configurations that present a reduction of shock wave strengths with respect 
to values produced in present configurations, and to validate results of 
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approximate analysis by means of experimental results. The investigation also 
is useful for developing satisfactory prediction methods. 

The investigation indicated that the analysis and experiments are in 

good agreement, and that sonic booms having the first shock strength of the 
2 

order of 1 lb/ft as predicted by the analysis are measured in the experi- 
ments. The investigation indicates some of the difficulties encountered in 
extrapolating near field data at large distances. The experiments have been 
performed at the Aeronautical Research Institute of Sweden, by the research 
group directed by Professor M. Landahl and Dr. G. Drougge, and have been 
performed by using a new experimental technique developed by the group ^ 

based on the higher order analysis developed by M. Landahl, I.L. Ryhming, 
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and others . The model has been designed by a research group at the Aero- 
space Laboratory of New York University, as a part of an investigation on 

3 4 

low boom configurations. ’ The latter group also performed the analysis of 
the experimental results described in this report. The work has been carried 
out under two NASA Grants: (1) NGR 52-120-001, assigned to the Aeronautical 
Research Institute of Sweden; and (2) Grant NGL 33-016-119, assigned to New 
York University. The experiments have been performed by H. Sorensen of the 
Aeronautical Institute of Sweden. The analyses have been performed by a 
team under the direction of Professors Ferri and Ting, which includes besides 
the authors, Mrs. F. Kung, Messrs. M. Sic lari and A. Agnone. This report 
discusses the results corresponding to an angle of attack close to the 
selected flight conditions for the SST, and at another slightly higher angle 
of attack. 


2 



II. Definition of the Problem 


Prediction of sonic boom signatures at large distances from the airplane 

5 6 7 8 

is based mainly on the Whitham second order analysis. * 5 5 The aircraft is 
represented by a distribution of multipoles and lifting elements derived from 
the airplane's geometry on the basis of linearized supersonic theory. The 
asymptotic solution is expressed in terms of the cross-sectional area distri- 
bution. S(x0) is a function of the distance x and azimuthal angle 0 of an 
equivalent body of revolution, and an equivalent cross-sectional area that 
represents an equivalent linear lift distribution. The pressure variations 
produced by a vehicle are strongly nonlinear near the vehicle; however, the 
nonlinear effects decrease rapidly moving away from the airplane. Therefore, 
at some finite distance from the airplane surface, linear theory is generally 
valid locally. Then the only nonlinear effects to be considered are the 
effects related to the variation of the speed of sound that when combined 
with the local variation of velocity of fluid produces distortions of the 
signal with respect to linear Mach waves. This second order effect is 
cumulative and therefore must be considered when propagation of distur- 
bances at a large distance is investigated. The Whitham theory takes into 
account these effects and permits us to determine the propagation of waves 
from a region where the disturbances produced by the vehicle are already 
sufficiently small, so that linear theory applies. 

The analysis permits us to extend the signal produced by the airplane to 
regions far from the airplane, provided that the disturbance produced by the 
airplane is known in a region where such disturbances are already small. 

The analytical determination of such disturbances is extremely difficult, 
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especially at high Mach numbers. The linear theory that is the most adapt- 
able type of analysis available for three-dimensional flow, cannot be used 
for the purpose of relating the shape of the airplane to the signature pro- 
duced in the near field, because it neglects higher order effects both in 
the definition of the local intensity of the signal, and in the decay and 
shape of the wave propagation. 

Experimental methods, where the wind tunnels are used as analog machines 

to define the signals at some distance from the vehicle, open a promising 

way to solve this problem because they do not have the limitations of the 

linearized theory. ^ Such approaches, however, have other limitations 

12 

that only recently have been recognized. The most direct approach for 
eliminating many of the uncertainties due to the simplifications introduced 
in the analysis to determine the actual signature of an airplane at some 
distance from the airplane, is to measure the complete signature experi- 
mentally. At Mach numbers on the order of 2 to 3, the disturbances produced 
by good L/D configurations at distances of the order of a few body lengths 
from the model are extremely small; therefore, linear theory applies here. 

9 

The experimental approach proposed by NASA , and used extensively for sonic 
boom determination, has been to build a model and to determine experimentally 
the disturbances produced at some distance from the airplane and then extra- 
polate the measured signal to the required height of flight by means of the 
Whitham theory. The concept is surely valid, and therefore presents one of 

the most direct ways of obtaining accurate information. However, it has been 
12 

shown that in applying the concept that more attention should be given to 
the type of measurements performed, and to the interpretation of the experi- 
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mental data, before such methods can give completely satisfactory results for 
configurations producing weak sonic boom signatures as required for future 
airplanes. One of the basic simplifications of the Whitham theory is that 
at large distances from the body, the details of the spanwise distribution of 
the sources of disturbances are not important provided that the span of the 
vehicle is very small with respect to the distance considered. Then, the 
airplane can be substituted by an equivalent axially symmetric body in any 
of the meridional planes around the axis of the body considered. This 
approximation has been introduced in the past, and used in the extrapolation 
of experimental measurements to large distances. Because of these simplifi- 
cations, the measurements have been performed only in the meridian plane 
normal to the plane of the airplane wing. These measurements have been used 
to determine by means of the Whitham theory, an axially symmetric body 
equivalent to the airplane configuration which is used to determine the sonic 
boom at large distances. This approach introduces two approximations: the 
first is that the region where the measurements are made is sufficiently far 
from the model, so that the disturbances are small; the second is that all 
dimensions normal to the axis of the airplane, including also the spanwise 
dimensions, are negligible with respect to such a distance so that the three- 
dimensional effects are neglected. Both assumptions are acceptable if the 
regions of the measurements are sufficiently far away from the model; however, 
the distance required for the validity of the first approximation is not the 
same distance required for the validity of the second. For airplanes having 
good performances, and wings of practical span, the first requirement can be 
satisfied at much smaller distances than the second. 
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In practice, distances from the model of the station where measurements 
are performed cannot be too large because of practical experimental diffi- 
culties. The precision of the measurements decreases sharply with an in- 
crease of distance because of the finite sensitivity of the instruments. 
Furthermore, the size of available wind tunnels and the size of the model 
required for the tests limits the ratio between distance of measurements and 
length of the model. Another important difficulty is due to the fact that 
the flow field produced by a wind tunnel is not absolutely uniform. In any 
wind tunnel, a nonuniformity of a few percent in Mach number exists. Such 
nonuniformities correspond to waves that interact with the wave pattern 
produced by the model. When these disturbances interact with the flow field, 
an error is introduced that is cumulative along the waves carrying distur- 
bances from the model to the plane of measurement. The strength of the 
waves that is produced by the model at a very large distance from the model is only 
slightly larger than the waves due to nonuniformities existing even in the 
most uniform wind tunnel. Therefore, the error introduced is proportionally 
larger at larger distances. This error tends to increase with Mach number. 

Thus, a compromise must be reached between accuracy required and distance of 

12 

the measurement where a single measurement is acceptable. 

The problem can be reduced substantially by the introduction of more 

complex techniques for obtaining experimental data, and extrapolating the 

results. Two different approaches have been proposed; one by the first 

12 

author of this report where the measurements are performed in a plane 
located at some distance from the model parallel to the plane where the 
sonic boom signal must be determined. The deviations of the stream surface 
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normal to the plane are measured in a region inside the shock generated by 
the front tip of the airplane. Such deviations are determined along several 
straight lines parallel to the flow direction to cover all of the flow in- 
side the shock. A stream surface is defined by the measurements that can be 
substituted for the airplane. Such a stream surface where the flow distur- 
bances are small is equivalent to the vehicle placed above. Then, the 

13 

Whitham analysis as applied by Walkden, is used directly to determine sonic 
boom at the required distance from this surface. The Whitham theory applies 
for this extrapolation provided that the disturbances are small, so that 
locally the linear theory is sufficiently accurate. In this type of an 
approach, the deviation of the stream surface is measured experimentally by 
measuring the stream deviation in a plane parallel to the ground in the 
entire region inside the front Mach cone. Therefore, the precision of this 
method depends on the compromise of two opposite requirements: (1) the 
accuracy of the linearized theory which increases with the distance from the 
body; and (2) the precision of measurements which decreases when the dis- 
tance increases. However, three-dimensional effects are accounted for 
accurately. The experiments presented here tend to indicate that a satis- 
factory compromise can be obtained for these two opposite requirements even 
if the measurements are performed at very small distances from the body. 

The present experiments have been performed by using the second method 
proposed by Landahl, Ryhming, Sorensen, and Drougge. 1 Here, the streamline 
deviation is measured for several streamlines starting on a cylindrical tube 
placed around the model having the axis parallel to the wind, and at small 
distances from the axis. In the experiments performed, the distance is 
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smaller than the length of the model. The deviation of each streamline of 
this tube is measured locally in several meridian planes. Two angles are 
measured: one gives the deviation in the meridian plane; and the second 
gives the deviation on the cylinder normal to the meridian plane. These 
data are used to determine at the axis of the cylinder the F function that 
is used in the Whitham theory, by means of higher order approximation that 
takes into account second order terms in the disturbance components and 
higher order terms in the curvature according to a theory developed by the 
proposers of the method. 

This method permits us to use smaller distances from the models than 
the other; however, it requires differentiation of the measured deviation 
which is difficult to do. In addition, the analysis assumes that the three- 
dimensional and thickness effects are small so that the disturbances can be 
extrapolated at the axis of the body. This last condition can produce F 
functions that are not singly valued. Both approaches are improvements with 
respect to the standard method, especially at Mach numbers of 2 or 3. 

III. Experimental Techniques 

The experiments were conducted in the Trisonic Tunnel of the Aeronau- 
tical Institute of Sweden, FFATUM 500, at Mach number 2.718.^ The tunnel 

2 

has a square test section of 50 x 50 cm with perforated walls for the 
transonic speed range and a flexible wall nozzle which allows the Mach 
number to be varied continuously between 1 and 4. It is a blowdown tunnel, 
which may be operated with a stagnation pressure up to 12 atmospheres, and 
a stagnation temperature range of 300°K - 400°K. A schematic design of the 
tunnel is shown in Fig. 1. 
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Pressure measurements were performed on the model at 2.6° and 3.2° 
incidence at two positions along the tunnel axis. The flow field measure- 
ments were conducted at two radial distances from the model axis correspond- 
ing to 0.271 and 0.558 times the length of the model in meridian planes 
spaced at 5° intervals from the plane of symmetry in the range between 0° 
and 90°. The meridian planes are defined by the angle 9 with respect to the 

plane of symmetry. The pressures were recorded almost simultaneously, since 

-4 

the time between the individual measurements were 1-10 sec. Schlieren 
photographs were taken of the flow field generated by the model and the 
pressure probe. 

The absolute level of accuracy of the results is very difficult to 
establish, because of the combined effects of the many possible sources of 
error. A number of precautions were taken, however, to reduce the magnitude 
and probability of significant errors. The facility instrumentation con- 
sists primarily of higher sensitivity pressure measurement devices for de- 
termining both stagnation and reference pressures. These pressures were 
calibrated carefully preceding the investigation. The free-stream pro- 
perties are considered accurate within the following limits: 

i 0.01 

p + 0 . YL 

t , 00 — 

The precision with which local flow quantities can be determined is es- 
timated by the Laboratory to be as follows : 
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Errors at M =3.0 


M 


1 


+ 0.07 



0 

<J 


= < 
= < 


+ 1.0 % 
+ 0 . 10 ° 
+ 0 . 10 ° 


The values of the errors in angles quoted here do not include the influence 
of the nonuniform flow on the probe. As it will be discussed later, the 
interaction of the shock with the subsonic flow in front of the probe, 
produces locally large errors; therefore, there such a type of measurement is not 
accurate. In addition some small influence due to Mach number gradients has 
been found experimentally. A report on such an influence will be made 
available later by the FFA group. However such effects are small and do not 
affect the results presented here. In order to analyze such effects, a 
comparison between measured and calculated values was performed for the front 
part of the model that is axially symmetric. The calculations performed by 
means of a numerical program based on a three-dimensional characteristic 
analysis is as accurate as the numerical approximation. The comparison 
that will be presented later, indicates that the estimated value of the 
error for g and a is larger than the maximum differences found between the 
measured and calculated values. 

The definition of the angles q and e, the support of the model, and the 
dimensions and point for rotation of the model are shown schematically in 
Fig. 2. The angles <j and q were measured by means of a probe as shown in 
Fig. 3. 

IV. Description of the Model 

In references 3 and 4, configurations have been investigated that 
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according to approximate analyses should produce sonic booms much lower than 
present configurations for equivalent conditions of flight. In order to 
determine the validity of such conclusions, and in parallel to experiments 
with the more accurate experimental techniques, one of these configurations 
has been selected for the tests. The configuration selected is shown in 
Fig. 4, where the dimensions correspond to the model tested. The wing is 
swept back at 72°. The wing profile has 2% thickness and is a symmetrical 
circular arc profile. The fuselage shape has a circular cross section; 
detailed dimensions of the fuselage area as a function of the distance are 
given in Table I. 

In the analysis it has been assumed that the fuselage is sharp at the 
front and rear tips, and that the wing profile is also sharp at the trailing 
and leading edges. The construction of the model had required some modifi- 
cation on the wing leading edge and fuselage front tip, and on the rear part 
of the fuselage. The modification introduced at the leading edge is required 
in order to avoid local separation. The modification at the rear part of the 
fuselage is required because of the presence of the support. The support 
increases the equivalent area in the rear part of the vehicle, and there- 
fore affects the sonic boom because it changes the final value of the equi- 
valent area distribution. 

In order to eliminate such effects, the equivalence between lift and 
cross-sectional area has been utilized, and a correction on the planform of 
the wing has been introduced. The area of the wing has been reduced in the 
region where the fuselage cross section is different from the design. The 
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redaction of wing area has been calculated to balance the increase of cross- 

sectional area from the reduction of lift and volume for a C of the wing 

L 

of 0.060, at M « 2.70. 

The design of the model is shown in Fig. 5, where the important differ- 
ences between the model used in the analysis and in the tests are shown. 

The configuration does not have a vertical tail because of the advisability 
of constructing a simple model. The volume due to the vertical tail has been 
included as an equivalent increase in volume of the fuselage. Actually, the 
vertical tail increases the equivalent length and is useful in decreasing 
the intensity of the rear shock. This effect is not important for this set 
of tests because the main interest is directed to the determination of the 
front shock, in view of the difficulty of obtaining representative measure- 
ments from the rear shock. 

The configuration selected does not represent a realistic airplane 
configuration. However, the total equivalent area distribution produced by 
the model tested at the selected C , has been selected from a more realistic 

Xj 

airplane configuration where the requirements of the volume and lift distri- 
bution have been selected from practical configurations. This configuration 
is shown in Fig. 5b. There are two reasons for the selection of the different 
configurations shown in Fig. 5, corresponding to the same total equivalent 
area. The first is of an analytical nature. The configuration shown in 
Fig. 5b is more complex* therefore, an approximate analysis that takes into 
account three-dimensional and nonlinear effects is not possible. Then an 
iteration procedure would be required, where the wind tunnel is used as an 
analog machine to select accurately the details of the configuration that 
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correspond to the equivalent area selected. The second is of an experimental 
nature; a minimum fuselage diameter was imposed by the balance. This required 
a reduction of the wing thickness and in increase of fuselage volume. 

Once a configuration has been obtained that gives low sonic boom, criteria 
are available to introduce changes that permit selecting an equivalent practical 
configuration similar to the configuration shown in Fig. 5b. 

V. Possibility of Obtaining Information on the Rear Shock from Wind Tunnel 
Tests 

The measurements in a wind tunnel of the rear shock produced by a model 
of this type are not indicative of the actual phenomena because of the lack 
of similarity. The presence of the sting, the difference in the size of the 
wake, and the lack of representation of the engine jets introduce severe 
errors in the intensity of the rear shock. The engines produce jets that mix 
with the outside flow producing a displacement thickness equivalent to sub- 
stantial area changes. Figure 6 indicates the equivalent area variation due 
to the jet downstream of one of the engines for an SST configuration which has 
four General Electric engines fully expanded. The total equivalent cross- 

2 

sectional area for a complete airplane at cruise is of the order of 800 ft . 

The difference between flight conditions and tunnel conditions in wake dis- 
placement thickness produces effects of the same order. These effects can 
be included in the analysis and have been included in the results presented 
in reference 4; however, they are difficult to simulate in small scale ex- 
periments of the type described here, and require a more complex research 
program than is required for the measurements of the front shock. 
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VI. Results of the Analysis 


The configuration presented has been analyzed for flight conditions 

at cruise for an airplane having a length of 300 ft in full scale. An 

2 

airplane of this dimension has 10,000 ft of wing area. The volume of the 

3 

fuselage is 44,000 ft , while the wing span is 58.8 ft and the wing 
3 

volume 10,000 ft . 

The analysis of the sonic boom has been performed for the conditions 

corresponding at cruise at 60,000 ft, M = 2.70, and total weight eaual to 

460,000 lbs. The analysis has been performed following the method outlined 
4 

by Carlson , which transforms the lift in an equivalent area. The lift has 
been determined by means of linear theory with some approximate second 
order corrections. The extrapolation of signature to ground level has been 
performed by means of the variable density program developed in reference 15. 
Figure 7 indicates the contribution of the equivalent area distribution of 
the different components. Table 1 gives the same information. The sonic 
boom signature obtained from the analysis is shown in Fig. 8 where four 
meridian planes are considered. The cut-off angle is around 60°. For 
comparison, the constant pressure analysis for average pressure equal to: 


^average V sea level ^60,000 ft 

has also been performed. The results are shown in Fig. 9 where the reflec- 
tion coefficient K is equal to 1.8. The analysis indicated that the front 
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shock for the weight assumed is of the order of 0.94. The ground reflection 
coefficient used in the analysis is 1.8. The area curve corresponding to 
the sonic boom in Fig. 8, is a smooth curve corresponding to sharp leading 
edges for the wing. The wing of the actual model produces a detached shock 
because the leading edge has been rounded. Therefore, some perturbations in 
the equivalent area distribution will be produced by this change. It is 
difficult to predict the actual intensity of these perturbations; however, 
the qualitative trend of the effects of such disturbances can be predicted. 
The round leading edge will produce initially stronger pressure rises than 
the sharp leading edge followed by an over expans ion. This is equivalent to 
a more rapid increase of the variation of the cross-sectional area in the 
region where the wing starts followed by a more gradual variation. The wing 
starts roughly at 50% of the length. The variations of the type of curves 
indicated in a and b of Fig. 10 can be expected. The sonic boom signature 
for these two distributions has been calculated. The results indicate that 
these local small perturbations tend to produce signatures that have two 
shocks of about the same intensity, Fig. 11. 

The calculations have been performed initially for a weight of 460,000 
lbs, at M = 2.70. The experiments have been performed at M = 2.72, and a 
C = 0.055, and C = 0.066 corresponding to the angle 2.6° and 3.2°. These 

ti 

conditions would correspond to a total airplane weight of 430,000 lbs and 
520,000 lbs. Then the condition of C = 0.055 is close to the design con- 

J_i 

ditions. The sonic boom for these two conditions has been calculated and 
is shown in Fig. 12. The effects of small local perturbations as shown 


15 



in Fig. 13 are also indicated. 

VII, Review of the Experimental Data 


The experimental data made available for the analysis are presented in 
Figs. 14 to 19. Figure 14 presents the measured values of 0 at r/L = 0.271 
for different values of 0 , while Fig. 15 presents the values of a for the 
conditions. Figures 16 and 17 present the same quantities for the distance 
r/L = 0,558. For several values of 0 , measurements are available for more 
than one position of the model along the axis of the tunnel. Figures 18 and 
19 present the results for 0=0 and r/L = 0.272 and 0.558 for the different 
positions. The figure indicates that the change of position does not affect 
the experimental results, giving an indication of the uniformity of the flow. 
Similar configurations on $ and e at the two distances but for a = 3.2 are 
given in Figs. 20, 21, 22, 23. As shown in the figures, the disturbances at 
— = 0.558 and a = 2.6° are extremely small; therefore, here linear theory 

J_i 

applies for this case. 

In addition, Schlieren pictures are available for all of these con- 
ditions. Figures 24a and 24b give the Schlieren photographs at ep = 0 and 
Cp = 90, for cc = 2.6°, and Figs. 25a and 25b for cl = 3.2°. The photographs 
permit us to locate accurately the position of the shocks, and therefore 
help in the interpretation of the experimental results. 

VIII. C orrections Introduced on the 0 and Justificat i on for the Corrections 
The probe that measures the angles e and <j has a diameter of 3 mm and 
has a spherical nose. The probe produces a shock and the flow region down- 
stream of the shock near the orifices is subsonic. The calibration of the 
probe is based on the assumption that the flow in front of the shock is 


16 



fairly uniform and the angle is determined on the basis of the assumption that 
a small inclination of the flow direction with the axis of the probe corres- 
ponds to the equal rotation of the pressure distribution around the axis of 
the probe. This assumption is correct everywhere if the flow is continuous 
because the probe is small. However, the situation changes when an external 
shock interferes with shock of the probe in the subsonic region. Consider 
Fig. 26 where the probe and the shock in the external flow are represented 
schematically. Interference with the calibration starts when the probe is 
behind the external shock, in position a, and continues until the orobe has 
crossed the shock and is at position b. 

Experiments have been performed**’ where this phenomenon has been in- 
vestigated at higher Mach number. In Fig. 27 taken from reference 16, the 
effects on the pressure distribution due to the interaction are shown. The 
variation of pressure distribution, due to the interaction when the sphere 
is in front of the external shock is very large (case c of Fig. 27). In 
this case if the pressure distribution is used to determine angles on the 
basis of a uniform flow calibration, then the measurement would indicate 
incorrectly that a large flow deviation exists in front of the probe. When 
the sphere is behind the external shock, then the interference produces an 
opposite effect, decreasing the value of the actual deviation. Because of 
these effects, the points measured in the vicinity of the shocks have been 
discounted and the diagram shown in Fig. 28 has been used for the analysis. 

In this diagram the position of the shocks have been determined very accu- 
rately from the Schlieren pictures. The peak deviation behind the shock 
has been obtained from the measurements of the shock inclination, and the 
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flow properties in front of the shock. 

In Figs, 29 and 30, the value of g calculated by means of three- 
dimensional characteristics are shown, and compared with the experimental 
data. The differences are extremely small and can be due either to small 
differences in geometry between the model and the shape used for the cal- 
culation or to measurements or calculation accuracy. Because the data of 
the sets of tests agree well, the first assumption seems to be more probable 
In order to determine the importance of some of the details of the 
pressure distribution, two alternate curves, 1 and 2 of Fig. 28 have been 
used in the analysis. Curve 1 considers only one shock, while curve 2 con- 
siders the existence of two shocks. The second shock is produced by the 
discontinuity in slope between fuselage and support. 

IX . D escription of the Method of Analysis Used 

The method used for the analysis has been described in reference 1 and 
has been developed by M. Landahl, I. Ryhming, H. Sorensen, and G. Drougge of 
the Aeronautical Research Institute of Sweden. In order to simplify the 
reading of this paper, the important features of the method are summarized 
here by reproducing a section of the report of reference 1. 

According to reference 2, the perturbation velocity components at large 
distances from a three-dimensional body in supersonic flow are given to 
second order by 


u 9 = — 

2 P 


u (x ,r ,e ) 
o’ o 5 o' 


( 1 ) 
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2 



v(x ,r ,e ) 


( 2 ) 
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w 2 = v(x o> r o ,9 o ) = 


(3) 


1 



where u^, v^» are the second-order components referred to a cylindrical 
coordinate system, 0 q is the meridian plane angle and r^ is the radial 
distance from the wind axis. 


x = x - Kr v - M 2 cp + Kr ^ 
o o ^ o se 


(4) 


r = r (1 - Ku - £ ) 

o P 30 
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and u, v, w are the corresponding values according to linearized theory and 
their potential is p. For large distances the following expansions were 
shown to hold: 
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where k = k 2 g 


The quantities measured are the flow deflection angles e and < 7 , along lin^s 
of constant r but in different azimuthal planes 0 = const. The angle e is 
related to v 2 as follows: 


e 


tan 


-1 



(13) 


or since u 2 = - v 2 /g to lowest order 


v 2 = (1 - ®- e + • • • 

From the aximuthal deflection angle $ we obtain 


(14a) 


w 2 w a 

or cp cp * ^ ra (14b) 

9 6 o 

Thus, cp can determined directly from the measurements, and by aid of 
®o 

numerical differentiation of a(0) one can also calculate cp « • From the 

o o 

measurements of q (x) , cp 2 can be determined by numerical integration, since 
to the order considered 


- - * r 


e (x' )dx' 


(15) 


r=const 


shock 


One can next obtain v from (2): 


v ~ (1 - M 2 u) (1- Ku) v 2 =- (1 + M 2 £■) (1 + K |) v 2 


(16) 
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Furthermore, (5) gives that 


r - r (1 - - e) (17) 

° P 

We have now the quantities needed to calculate F from (7). The Mach line 
parameter y and the angle parameter 0 q , finally, can then be determined from 
(11) and (12). 

For an axisymmetrical flow field w is zero as well as 0-derivatives. 

The measurements and the evaluation of the F-function is then simplified, 
since it is only necessary to consider a single azimuthal plane. 

X. Application of the Method to the Present Experiments 

The method presented is an improvement with respect to the methods 
used previously, because it takes into account nonlinear effects which are 
important in the near field and because it introduces corrections for the 
near field effects due to three-dimensional phenomena. In applying such a 
method to the present experiments, two different difficulties have been 
encountered. 

The first difficulty is related to the evaluation of cp in equation 7 

90 

of the analysis, which required differentiation with respect to 0 of the 
measured quantity q. Such an operation is difficult to perform accurately, 
because small errors on the measured values of a introduce variation of the 
value of cp . An attempt has been made first to perform a Fourier analysis 
of all the measured values, between 0=0 and 90°. However, for the front 
part the experimental accuracy is insufficient, while for the region of the 
wing where the contribution of cp. Q to the final results has significance, 
the series does not converge. Therefore, a different approach has been used 
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in the analysis. In the region of the wing the values of a between 0 and 
30° have been plotted as a function of 0, and a curve has been drawn connect- 
ing the experimental points. Then the value of at cp = 0 has been 

obtained from the curve. For the region where the model has axial symmetry, 
the values of So7de given by the analysis have been used. The contribution 

of cp to the F function in this region is negligible. Therefore such an 

09 

approach is justified. However the contribution of the term is important 
in the region of the wing indicating the importance of the three-dimensional 
effects . 

The second problem encountered in the determination of the F function 
is related to the existence of double values of the F function due either 
to the three-dimensional effects or to rapid expansions (corners) at the 
surface of the body.^ This effect is due to the fact that the F function 
is calculated at the axis and can be understood clearly if a simple configura- 
tion is analyzed. Consider, as an example, a cone cylinder at zero angle of 
attack. At the expansion corner, a family of expansion waves is produced. 

Fig. 31. If these expansion waves are extrapolated to the axis, then two 
sets of values corresponding to the same y for the F function in the region 
AB of the figure. Similar effects exist when the model is three-dimensional. 
Then waves generated from different points of the model at different r reach 
the same point in the meridian plane analyzed. In order to obtain some 
criterion for the solution of this difficulty, a detailed analysis has been 
performed for the case of the cone cylinder of Fig. 31 for the case, the 
actual stream deviation e at several values of r can be numerically obtained 
by means of the method of characteristics. Then it is possible to deter - 
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mine the F function that corresponds to the e distribution obtained at some 
distance from the axis r and compare with the F distribution obtained by in- 
troducing the discontinuity. Such analysis indicates that the double value 
should be eliminated by adding and subtracting equal areas in F on the 
Mach wave from the corner as discussed in reference 7, (see Fig. 31), in 
the region where multiple values exist. The F function obtained from the 
data is shown in Fig. 32. Figure 33 shows the F function used in the analy- 
sis. The same method of analysis has been used for the data at r/L = 0.271, 
and for a - 3.2°. The F function obtained from the data at 0.558 has been 
used to obtain the sonic boom corresponding to an airplane 300 ft long flying 

at 60,000 ft, for the C of 0.055 which corresponds to a total weight of 

L 

430,000 lbs. Figure 34 gives the sonic boom shape obtained from the F 
function by using the variable pressure program, and Fig. 35 the same re- 
sults obtained by using constant pressure program corresponding to the 
square root of the two values at the flight and ground level. The amplifi- 
cation due to the reflection of the ground is assumed to be equal to 1.8 
the incoming signal, and is included in the values presented. In order to 
determine the importance of the three-dimensional and nonlinear effects, 
the calculations have been repeated, by applying directly the Whitham theory 
to determine the F function, and neglecting all the higher order terms, 
curve a and only the three-dimensional effects (cp ) curve b. Figure 36 
shows the sonic boom obtained, when these quantities are neglected. The 
figure indicates that the three-dimensional effects are important at least 
when r/L is small. A similar calculation has been performed for the data at 
a = 2.6° and r/L = 0.271. The two e distributions for the analysis are as 
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shown in Fig. 37. The results obtained are shown in Fig. 38. They are very 
similar to the results of Fig. 35. The position of the second shock is 
different. 

The same analysis has been performed for the data corresponding to 
a - 3.2 e . This condition corresponds to a vehicle that has a length of 
300 ft, weight of 530,000 lbs, flying at an altitude of flight of 60,000 ft. 

The results are shown in Figs. 39, 40, and 41. Figure 39 gives the distribu- 
tion of e at r/L - 0.558 used in the analysis. Figure 40 gives the F function 
calculated and the curve used after being corrected by eliminating the double 
value, and Fig. 41 gives the shape of the sonic boom derived from the F function. 
XI. Conclusions 

The experimental investigation performed permits us to reach the following 
conclusions : 

2 

1. Sonic booms having peak values of the order of 1 lb/ft as 
predicted analytically in reference 3 have been measured. The sonic 
booms obtained have near field signatures as predicted in reference 3 
for cruise conditions and in reference 17 for the acceleration phase. 

The distribution of equivalent cross-sectional area tested corresponds 
to an airplane shape that has the volume, length, and lift requirements 
of a practical airplane configuration; however, substantial additional 
work is required to investigate if all other aerodynamic requirements 
related to a practical configuration can be met. 

2. The nonlinear and three-dimensional effects are of primary importance 
for the determination of the correct values of the sonic boom from 
measurements at small distances from the model. 

3. More complex experimental techniques where such effects are determined 
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are required when near field measurements are made as recommended in 
references 1 and 12. 

4. The experimental method proposed in reference 1 gives satisfactory 
results. 

5. Improvements are still required in the experimental techniques 
and in the analysis in order to measure and determine with better 
accuracy all of the required quantities. 
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Schematical view of test set-up 
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Fig. 3 Design of yaw probe 
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Fig. 5a Model design 
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Fig. 6 Variation of Total Equivalent Area Due to the Jet Exhaust of an Engine 
as a Function of Distance from the Exit of the Jet. 




Fig. 7 Equivalent area distribution 
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Fig. 9 Sonic boom signature constant pressure program*? K ^P 60)000 = 1017 lbs/ft 
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CURVE a OF FIG. 13, Q.* 0.066 
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Fig. 12 Sonic boom signature at M = 2.718 L = 300 ft, h = 60,000 ft, = 430,000 lb (t = 0.055) 
W„ = 516,000 lb, ( <t = 0.066) 
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Fig. 14d Continued 
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48 Fig. 14e Continued 
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Fig. 15a Experimental values of a as function of distance at several meridian planes 
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Figo 15b Continued 
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Fig, 16b Continued 
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Fig. 16d Continued 
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Fig. 17b Continued 
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Fig. 20c Continued 
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Fig. 20f Continued 
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Fig. 22e Continued 
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Fig. 23a Experimental values of a as function of distance at several meridian 
planes , a - 3.2°, r/l Q = 0.558 
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Fig. 25b Schlieren photographs a 





Fig. 27 Pressure distribution on a hemisphere, M = 4.6 interacting 
with a shock that deflects the flow of 5 . Arrow indicates 
point at which disturbance meets model surface. 
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Fig. 29 Comparison between experimental and calculated values 
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Fig. 31 Double values of F (y) for cone cylinder 
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Fig. 34 Sonic boom signature at r/L =» 200, for L = 300 ft and L = 0.555, variable density program 
(Curve 2 of Fig. 28) 0 ° ° 
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Fig. 35 Sonic Boom Signature (constant pressure) at r/L - 200 Kp = 1017 lb/ft 
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Fig. 36 Comparison of Sonic Boom Signatures for several types of analyses r/L = 200 
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Fig. 37 Measurements of deviation angle of SST model at r/L Q - 0.271, ( t= 0.055, a _ 2.6 at several 


X(mm) 


positions 




Fig. 38 Sonic boom signature for a = 2.6 


C = 0.055 at r/L = 
o 


200 from data at r/L = 0.271, Kp - 1017 lb/ft 
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Fig. 40 F(y) curve for ot = 3.2° from data at r/L Q = 0.558 
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Fig. 41 Sonic boom signature for a = 3.2°, ^ = 



